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ABSTRACT 



cd It has been widely assumed that partially quenched chiral perturbation the- 

I ory is the correct low-energy effective theory for partially quenched QCD. Here 

^ we present arguments supporting this assumption. First, we show that, for 

^ partially quenched QCD with staggered quarks, a transfer matrix can be con- 
structed. This transfer matrix is not Hermitian, but it is bounded, and it can be 

^ used to construct correlation functions in the usual way. Combining these obser- 

QQ vations with an extension of the Vafa-Witten theorem to the partially quenched 

■^ theory allows us to argue that the partially quenched theory satisfies the cluster 

^ property. By extending Leutwyler's analysis of the unquenched case to the par- 

• tially quenched theory, we then conclude that the existence and properties of the 

(^ transfer matrix as well as clustering are sufficient for partially quenched chiral 

C^ perturbation theory to be the correct low-energy theory for partially quenched 

Tl QCD. 
> 

X 
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I. INTRODUCTION 

Partially quenched chiral perturbation theory (PQChPT) has been extensively used in 
the analysis of numerical computations of hadronic quantities in lattice QCD. In such com- 
putations, one has the freedom to vary valence quark masses (masses of quark operators 
appearing explicitly in correlation functions) and sea quark masses (masses of quarks ap- 
pearing in the fermion determinant of the theory) independently. This generalized version 
of QCD, which is commonly referred to as partially quenched QCD (PQQCD), contains full 
QCD as the special case in which valence and sea quark masses are set equal to each other 
(for each flavor) [1]. PQChPT is, correspondingly, the generalization of chiral perturbation 
theory to the partially quenched setting.^ 

The ability to vary valence and sea quark masses independently is useful for a variety of 
reasons. First, the computation of quark propagators needed for the contractions making up 
a correlation function is significantly less expensive in most applications than the generation 
of gauge field configurations, which depend on the sea quark mass. With fixed computational 
resources, it can thus be an advantage to generate data for a number of valence quark masses 
on an ensemble of gauge configurations with a given sea quark mass. 

Second, PQQCD contains full QCD with the same set of sea quarks. It follows that the 
low-energy constants (LECs) of the effective theory for the partially quenched theory are 
those of the real world, because, by definition, the LECs do not depend on the quark masses 
[3].^ It turns out that in a number of cases, it is easier to determine these LECs by varying 
the valence and sea quark masses independently; having more parameters to vary provides 
more "handles" on the theory. 

Third, on the lattice, PQQCD can be generalized to QCD with a "mixed action," in which 
not only valence and sea quark masses are independently chosen, but also the discretization 
of the Dirac operator is different for valence and sea quarks [4] . The continuum limit of such 
a theory is a partially quenched theory; a fine tuning is generally required to make valence 
and sea quark masses equal in the continuum limit. 

Fourth, lattice theories with staggered fermions that use the fourth-root procedure [5] 
to ehminate unwanted "taste" degrees of freedom are closely related to PQQCD because 
there is a mismatch between the (rooted) sea quarks and the (unrooted) valence quarks 
[6-8]. PQChPT plays an important role in the development of the effective chiral theory for 
rooted staggered fermions [9, 10], which in turn provides evidence that the rooting procedure 
is indeed valid. For more discussion of these issues see Refs. [7, 11-13] and references therein. 

However, it is less clear than in the case of full QCD that the partially quenched chiral 
theory is indeed the proper low-energy effective field theory. This is because PQQCD violates 
some of the properties of a healthy quantum field theory: The path integral definition of 
PQQCD includes an integral over ghost quarks, which have the same quantum numbers 
and masses as the valence quarks, but which have bosonic, rather than fermionic, statistics. 
The partially quenched theory thus violates the spin-statistics theorem. The reason for the 
presence of ghost quarks is that their determinant cancels the valence quark determinant 
[14]; it is this very cancellation that makes the gauge configurations independent of the 
valence mass. 

In Ref. [15], Weinberg conjectured that the validity of chiral perturbation theory (ChPT) 



^ For a review, see Ref. [2] . 
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as a low-energy effective tfieory for tlie Goldstone sector of QCD follows from the basic 
properties of a healthy quantum field theory, which include analyticity, unitarity, cluster 
decomposition, and symmetry considerations. It was assumed that the S'-matrix calculated 
with the most general local Lagrangian consistent with a certain symmetry group is the 
most general possible S'-matrix consistent with these basic properties. This was then used 
as a starting point for the systematic development of ChPT to obtain S-matrix elements 
as an expansion in terms of the pion momenta and masses, following a well-defined power- 
counting scheme. The implicit reliance of this argument on unitarity, though, appears to be 
an obstruction to extending this line of reasoning to the partially quenched case, which is 
certainly not unitary. 

A justification for ChPT as the low-energy effective theory for QCD based on a somewhat 
different set of arguments was presented by Leutwyler [16]. In this justification, the most 
important ingredients, in addition to symmetries, are locality and the cluster property of 
the underlying theory (full QCD), while unitarity is not used. Locality and clustering guar- 
antee the existence of vertices in the effective theory that are independent of the correlation 
functions in which they appear, and, consequently, the existence of a loop expansion. We 
have found this approach more useful for the case of PQQCD than that of Ref. [15]. By 
construction, PQQCD is local. The question then becomes whether the theory also satis- 
fies the cluster property. The main goal of the present article is to collect the theoretical 
evidence that this is indeed the case. Our argument will be based on three ingredients: 
the existence and properties of a transfer matrix for PQQCD, the observation that chiral 
symmetry breaking is expected to take place in PQQCD because it takes place in the full 
theory contained within PQQCD [17], and an extension of the Vafa-Witten theorem [18] 
about the absence of spontaneous symmetry breaking of vector-like global symmetries (this 
part of the argument extends the argument already given in Ref. [17]). The combination 
of these ingredients allows us to argue that the partially quenched theory satisfies the clus- 
ter property, under a set of mild additional assumptions similar to those used in Ref. [16]. 
With this result in hand, the justification for PQChPT as the low-energy effective theory 
for PQQCD then follows, much like it does for the case of full QCD. For technical reasons, 
we limit ourselves to lattice QCD with staggered quarks, but we believe that the extension 
to other discretizations of QCD is relatively straightforward.'^ An earlier account of part of 
this work appeared in Ref. [20]. 

This article is organized as follows. In Sec. II we construct the transfer matrix for PQQCD 
with staggered quarks. We find that the transfer matrix is not Hermitian, but is nevertheless 
bounded. It turns out to be instructive to consider the free theory in some detail, and this is 
done in Sec. III. In particular, the free theory clearly demonstrates that unitarity is violated 
in PQQCD. 

Then, in Sec. IV, we turn to the effective theory, PQChPT. In Sec. IV A, we argue 
that the extension of the Vafa-Witten theorem about vector-like global symmetries in QCD 
implies (up to certain mild assumptions) that the cluster property holds in PQQCD. We 
discuss the role of rotational symmetry in Sec. IV B, and synthesize all our observations 
into an argument for the correctness of PQChPT as the low-energy effective theory for 
PQQCD in Sec. IV C. For technical reasons, we need to assume that the pion masses remain 
real in PQQCD, despite the fact that the corresponding transfer-matrix Hamiltonian is not 
Hermitian. This assumption is strongly supported by numerical evidence. In Sec. IV D, 
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we use CPT symmetry to argue that all low energy constants in PQCliPT have the same 
phases (real, with usual conventions) as they would in a Hermitian theory. This supports 
our assumption that the masses are real. Nevertheless, as we show in Sec. IV E, in the 
nonstandard case where ghost-quark masses are not degenerate with valence-quark masses, 
the effective theory shows that complex masses may arise in some ranges of quark and ghost 
masses, as may a new phase transition. This is mainly of academic interest, because only the 
case of pair-wise degenerate ghost and valence-quark masses is relevant in describing lattice 
gauge simulations. Our conclusions are contained in Sec. V. There are two appendices, both 
concerned with the free theory. In App. A we show completeness of a basis of (right or left) 
eigenstates, and App. B discusses a path integral formulation of the free theory. 



II. TRANSFER MATRIX 

In this section, we construct the transfer matrix for a gauge theory coupled to fermionic 
and bosonic or "ghost" staggered quarks. Each of these quarks can have an arbitrary mass, 
but we will require all masses to be positive.^ In Sec. II A we construct the transfer matrix 
for ghost quarks in a background gauge field. In Sec. II B we discuss some properties of 
the ghost transfer matrix; in particular, we show that it is bounded. Then, in Sec. II C 
we combine this with the transfer matrix for a gauge theory with only fermionic quarks to 
arrive at the complete transfer matrix, and show that it is invariant under PT and CPT 
symmetry. We will use the "double time-slice" construction for both fermionic and ghost 
quarks [22, 23]. 

A. Ghost sector 

The staggered action is 

2 ^ Vt^ix) {xHx)U^,ix)xix + /i) - x"^(x + /i)[/^(x)x(a;)) + mx^ix)xix) 

(2.1) 



s=i: 



in which 

r/,(x) = (-l)-i+--.- , (2.2) 

and where x ci-nd x^ cire the staggered fields, Ufj_{x) are the link variables, and color indices 
are suppressed. We will denote a lattice gauge field consisting of all link variables f/^(a:) by 
U. If X and x^ are Grassmann, they are independent of each other, and x^ can then also be 
denoted as x, as is often done. However, here we are interested in the ghost-quark sector of 
the theory, for which we take x{x) to be a c-number. In order that the path integral over 
the ghost fields be convergent, we need to take xK^) ^o be the Hermitian conjugate of x{x). 
With this choice, the partition function 

Z(P() = JU dxldx. exp (-5) (2.3) 



The phase of a ghost-quark mass cannot be changed by chiral transformations [19, 21]. 



is well-defined, as long as we take m > 0, which we will assume throughout this article. 

In order to construct the transfer matrix representation of Z, we find it convenient to 
introduce real fields 0i and 02 through 

X{x) = r]i{x)(l)i{x) + #2(3;) , xH^) = mix)(l)i{x) - #2(2;) . (2.4) 

Choosing temporal gauge, U^i^x) = 1, and splitting x — )■ (x, t = X4), we rewrite S as^ 

S = ^<^i(^0i(x,t)02(x,t + l)-02(:r,t)0i(x,t + l)) (2.5) 

+i'^r]'j{x) ((f)i{x,t) ReUj{x,t)(f)2{x + j,t) + (f)2{x,t) Re Uj{x,t)(f)i{x + j,t)] 

.:E.(^l(-.(^.«..-,(^.0.(-I.O^..(^V,:..,(....(..J.,; 

j 

+m(0i(a;)2 + 02(a:)^) 

in which 

V'j{x) = rij{x)v4{x) , (2.6) 

and 

ReU,{x) = l(Uk{x) + U;ix)) , (2.7) 

ImUkix) = Y^{Uk{x)-U;{x)) . 

Next, we divide the lattice into even and odd time slices, and rename the fields $, respectively 
n, defining 

t = 2k: 0i(x,t) = <l>i,fc(x) , 02(x,t) = -$2,fe(x) , (2.8) 

t = 2A; + 1 : 0i(x,t) = U2,k{x) , 02 (^,t) = Ili,k{x) . 

Accordingly, the action can be rewritten as 



s = $^|x^^($i,.(x)ni 

k \ X 



h,k{^) + ^2A^)^2A^)) (2-9) 

+?/_[$i,fc, $2,fc;W(2A;)] + 'Hoi'^i, $2; m] 

+ J2\J2~^ (ni,fc(f)<l>i,fe+i(x) + Il2,kix)^2,k+l{x)) 
k K X 

+n4Ui^k, n2,fc; u{2k + 1)] + no[Ui, U2; m] 

^ The fields (fii and 02 are always transposed when they appear as the first factor in a bilinear in these 
fields. Wc also use that 774 (x+j) = — 774(0;). 



exp 



in which U{t) denotes the gauge field at a time shce t, and 

•H±[^i,^2; W(t)] = 5^|±5^zr/;(f)(^i(f)Re[/,(f,t)^2(x + j) + (1^2)^ 

X K j 

-^^iVjix) ('^i(x) Imf/j(f,t)^i(f + j) -(1^2) 
j 
?^o[^i,^2;m] = mJ2{'^i,k{^f + '^2,k{xf) . (2.10) 

X 

After shifting A; to A; + 1 in the H- term, we define a kernel 

Ti^i,k+i<^2,k+i; $i,fc$2,fc) = / n ^ni,fc(^) y n di^^^y) >< (2.ii) 

y y 

^i($i,fc(f)ni,fc(f)+<i>2,fc(f)n2,fc(f)-ni,fc(f)<i>i,fc+i(f)-n2,fc(x)$2,fc+i(:r)) 

+H- [$1,^+1, $2,^+1; W(2(fc + l)]+H+[n2,fc,ni,fc; U{2k + l)] 

+Ho[$i,$2;m]+Ho[ni,n2;m] 

The claim is then that 

T($l,fc+l$2,fc+i; $l,A:$2,fc) = ($l,fc+l$2,fc+l|TG,fc(W)|$i,fc$2,fc) , (2.12) 

with 

T (U) = e-^-[*i'*2'"(2(fc+i)]--Ho[<i>i,-J'2;m] g--H+[n2,ni;W(2fc+i)]-Wo[ni,n2;m] ^2 ]^3^ 

in which the Hermitian operators $a(a?) and na(x) obey the commutation rules 

[^a{x),flb{y)]=tSs,y-6ab. (2.14) 

This is proven by inserting a complete set of states into the right-hand side of Eq. (2.12): 
($i,fc+i$2,fc+i|TG,fe(W)|<l>i,fc$2,fc) = (2.15) 

X (ni,fcn2,fc|$i,fc$2,fc) 

= f dllik [ dU2k e-'^-[*i''=+i'*2.'=+i'"(2(fc+i))]g-w+[n2,fc,ni,fc;W(2fe+i)] 

X g*(*i,fc+ini,fc+'i'2,fe+in2,fc-<J'i,fcni,fc-<j>2,fen2,fc) 

where we omitted the explicit arguments x from the fields. Restoring these, the last line of 
Eq. (2.15) coincides precisely with Eq. (2.11). The ghost partition function in a gauge-field 
background U is now given by 



T/2 



ZG(W) = Tr I J]fG,fe(W) 
if T is the time extent of the lattice. 



(2.16) 



, fc=i 



B. Properties of the ghost transfer matrix 

The transfer matrix (2.13) is not Hermitian, but it is bounded. The proof is as foUows. 
The operator 'H_[$i, $2; W(2(/c + l))] +'Ho[$i, $2; itt] consists of a positive semi-definite part 
(if m > 0), and an anti-Hermitian part containing the gauge field. Moreover, it commutes 
with its Hermitian conjugate, because all ^a{x) commute among themselves. Therefore, the 
operator 

f iU) = e-^-[*i'*2;W(2(A:+l)]-Wo['J'i,*2;m] ^2 i^j) 

is normal and bounded,^ 

II TiiU) II < 1 . (2.18) 

A similar argument applies to 

T2(W) = (,-'H+in2,ni;U{2k+i)]-no[Ui,n2;m] ^ /2.19) 

and thus it follows that Tq itself is bounded, 

II faiU) II < II fi(W) nil f2{U) II < 1 . (2.20) 

This establishes that all eigenvalues of Tq have an absolute value less than or equal to one. 
If, moreover, the eigenvalue Aq with maximal absolute value is unique, correlation functions 
in this theory decay exponentially with distance. 

The transfer matrix (2.13) may be assumed to have a complete set of right and left 
eigenstates. This is equivalent to saying that, if we block-diagonahze Tq (put it in Jordan 
normal form), there are no blocks of the form 



A K 
A 



(2.21) 



with K 7^ (or generalizations of this form with higher degeneracies). We will see in Sec. Ill 
(and App. A) that this situation does not occur in the free theory. The condition that it 
happens on some nontrivial gauge field in the interacting case then restricts such fields to 
a subspace of co-dimension one (or more) in the full space of gauge-field configurations. 
Therefore, for "most" (meaning all gauge fields except a set of zero measure in the space of 
all gauge fields), the transfer matrix Tc,kiM) can be completely diagonalized, and complete 
sets of right and left eigenstates exist. 

There is a caveat, however. In the next section, we will incorporate the ghost transfer 
matrix we have constructed thus far into a transfer matrix for the entire theory, including 
quantized quarks and gauge fields as well as ghosts. In order to conclude that the total 
transfer matrix does not have any nontrivial Jordan blocks of the form (2.21) (or general- 
izations thereof, c/. App. A), we would have to construct the full Hilbert space for the entire 
matrix, something we do not know how to do.^ One might, however, consider a hybrid 



^ We use the Euclidean norm, which, for a matrix A, is defined as the positive square root of the largest 

eigenvalue oi A^ A. 
^ For the free theory it is obvious that it can be completely diagonalized, because in that case the Hilbert 

space is the direct product of the free Hilbert spaces for the quark and ghost sectors. 



construction of the entire theory, in which ferniions and ghosts are treated in a canonical 
formahsm, and gauge fields are taken into account through the path integral. What this 
means is that correlation functions with quarks or ghosts on the external lines are first con- 
structed in the transfer-matrix formalism, in an arbitrary fixed gauge-field background. The 
full QCD correlation functions are then obtained by integrating these correlation functions 
over the gauge fields. We first exclude from this integral the measure-zero set of gauge fields 
for which the Jordan normal form of the transfer matrix in a fixed gauge-field background 
may be nondiagonal. Since we then have a complete set of eigenstates in each background 
field, we think it reasonable to assume that the entire transfer matrix has a complete set of 
eigenstates. 



C. Full transfer matrix 

It is straightforward to combine the transfer matrix for ghost quarks constructed in 
Sec. 11 A with the transfer matrix for lattice QCD with staggered fermions. First, from 
Refs. [22, 23], the fermionic transfer matrix for a staggered quark in a fixed background 
gauge field lA can be written in the form 

T^^,(W) = e^'["(2('=+i))] e^t'"] e^m2k+i)] ^ (2.22) 

with B Hermitian. This translates the system from double time slice k {cf. Eq. (2.8)), with 
gauge fields U{2k) and U{2k + 1), to the next double time slice A; + 1, with gauge fields 
lA{2{k + 1)) and lA{2{k -|- 1) + 1). In more detail, the factor e^ takes care of the hop within 
the double slice, connecting slice 2k to slice 2^-1-1 (and it contains spatial terms on slice 
2/c + 1), whereas the factor e^ hops from slice 2/c + 1 to slice 2{k + 1) in the next double 
slice (and it contains spatial terms on the slice 2(fc + 1)). 

The ghost transfer matrix of Eq. (2.13) can be written similarly as: 

fa^iU) = e-'^-[*i'*2;W(2{fe+i))] ^-mm] ^-H+[Ui,n2;U{2k+i)] ^ (^2.23) 

where 

^-n[m] ^ g--Ho[l>i,<i'2;m] g--Ho[ni,n2;m] ^2 24) 

We can then write the transfer matrix for the total theory of QCD with staggered quarks 
and ghost quarks in the form 

f,„,,, = t^/2 e^' e-^-[*^'*^^"] fu e-^H e^ e^ e-^+P^'"^^"] f'J^ , (2.25) 

with Tu the transfer matrix of the pure gauge theory constructed in Ref. [24]. The gauge 
field transfer matrix hops between single time slices, and therefore a factor Tu needs to be 
inserted between e^e-^+[ni,n2;W] ^^^^ ^a g-w_[$i,<i>2;W]^ with matrix elements 

{U{2{k + l))\fu\U{2k + l)) . (2.26) 

Then, when one writes the partition function as a trace over a power of the transfer matrix, 
the factors Tjj combine to hop the gauge field from slice 2k to 2k + 1 (on the right of 
Eq. (2.25)), or from slice 2{k + 1) to slice 2{k + 1) + 1 (on the left of Eq. (2.25)). 



vr^iix) {vry = 


<^2(-x) ., 


VT^2{x){VTy = 


^li-x) .. 


vrtii{x){vry = 


-M-x 


vrn2ix) (vry = 


-ni(-x 



Although Ttotai is not Hermitian, and therefore is not required to have only real eigenval- 
ues, there are significant restrictions on the eigenvalues coming from discrete symmetries. 
Most importantly, Ttotai is invariant under the antiunitary symmetry PT, the product of 
parity and time-reversal symmetries. Under PT, the ghost fields in Eq. (2.25) transform 
according according to 



(2.27) 



where VT is the antiunitary operator that generates the symmetry. Note that, because 
it is antiunitary, PT it leaves the commutation rules, Eq. (2.14), unchanged. In the free 
case {Uj{x,t) = 1), the invariance of the ghost part of the transfer matrix. To, under PT 
can be easily checked using the definitions of 'H± and Tio, Eq. (2.10), and the relations 
rij{x + j) = rij{x) and tj'Ax + j) = —rj'Ax). In the interacting case, the ghost transfer matrix 

TciU) is of course not invariant under Eq. (2.27) for fixed background fields Uj{x, t); we must 
also let the time-slice gauge field operators transform. The transformation rule is standard: 

vr Uj{x) {vry = u^ji-x) = u]{~x - J) . (2.28) 

Of course the quark and pure gauge parts of Ttotai are also invariant under PT. 
Since Ttotai is invariant under PT, 

vr Ttotai {Vry = Ttotai , (2.29) 

each eigenvalue of Ttotai niust be either real or one of a complex conjugate pair: 

ftotall"^) = A|^) => ftotaliVri^)) = X*{Vr\^)) . (2.30) 

Non-Hermitian Hamiltonians with PT symmetry have been studied extensively by Bender 
and others [25]. PT-symmetric theories may often be redefined to give an acceptable unitary 
theory. Here, however, we do not want to make any such redefinitions, since the path 
integral is given, and the unitarity violations due to a non-Hermitian transfer matrix (or 
Hamiltonian) are as expected for a theory with spin-1/2 bosons. 

For future reference, we note that Ttotai is also invariant under charge conjugation sym- 
metry, C, and hence under the combined symmetry CPT, with 

cvr^iix) (cvry = <l'i(-f) , 

CVr ^2ix) iCVT)^ = ^2{-x) , 

CVTili{x) (CVry = -fli(-f) , (2.31) 

CVTMx) (CVT)^ = -M-x) . 

The gauge field transforms under CPT as 

cvr Uj{x) {cvry = uj{-x - J) , (2.32) 

with T the matrix transpose. 



From Eq. (2.31), we can find the CPT transformation rules for the c-number fields x{^) 
and x^{^)y which are, from Eqs. (2.8) and (2.4), the eigenvalues of (linear combinations 
of) the operators $i(a;), $2(x), tli{x), and 112 (a^) on each time slice. Note that, since the 
time-translation operator in Euclidean space is exp^—Hxi), we should not send x^ — )■ — X4 
under this symmetry if we want the Euclidean action, as opposed to merely the Hamiltonian, 
to be invariant.^ We then have 

CPT: xix,x,) ^ (-iy^x^^(-x,x,), 
CPT: xH^,x^) ^ i-iy^x^i-x^x^) , (2.33) 

where the factors of (—1)^"* arise from the minus signs in the last two equations in Eq. (2.31). 
It is straightforward to check that the action, Eq. (2.1), is unchanged by this transformation. 
For staggered quarks, the action is 






S = y2{oY] Vfi{x) {q{x)U^{x)q{x + /i) - q{x + /i)t/^(x)g(x)) + mq{x)q{x) \ , (2.34) 



where now q{x) and q{x) are independent Grassman-valued fields. Starting from the naive- 
quark transfer matrix Hamiltonian [23] or the "reduced-staggered" transfer matrix in [22], 
one can derive the transformation rules for q{x) and q{x) that correspond to Eq. (2.33): 

CPT: q{x,X4) ^ {-ly^ f{-x,Xi) , 

CPT: g(x,X4) ^ -(-l)"-'^ g^(-x,X4) . (2.35) 

The extra minus sign in the second equation in Eq. (2.35) (as compared to Eq. (2.33)) makes 
up for the minus sign coming from Fermi statistics when taking the transpose of the action. 
We end this section with a few comments. First, implicitly, we have only considered one 
flavor of quarks and one flavor of ghost quarks. The generalization to arbitrary numbers 
of each is immediate. In addition, we have generalized beyond PQQCD by choosing an 
arbitrary (positive) mass for each quark or ghost quark. In PQQCD, each ghost quark mass 
is equal to the mass of one of the fermionic quarks, thus turning that quark into a valence 
quark. Sea quarks appear as fermionic quarks without ghost partners. A further issue arises 
from the use of the staggered action (without rooting) for the ghosts and quarks, which 
implies that each flavor comes in four tastes. This is actually not a serious restriction for 
PQQCD, since extra (unwanted) species of valence quarks and ghosts are harmless: they 
have no effect on any processes if we choose not to put them on external lines (in the standard 
case where valence quarks and ghosts are degenerate). Further, if one wishes to avoid extra 
tastes in the sea, it would be completely straightforward to use any alternative discretization 
for the sea quarks that has a transfer matrix, such as unimproved Wilson quarks. 

III. THE FREE THEORY 

This section focuses on the ghost transfer matrix, Eq. (2.13), in the free theory, i.e., in 
the case where the background gauge field Uj{x,t) = 1. We work in the limit of vanishing 

® There is however a linear time-inversion symmetry of the action in Euchdean space, which is simply the 
time-direction equivalent of the spatial inversions. We focus instead on the antiunitary symmetry of the 
Hamiltonian because it will be useful to us in constraining the chiral theory. 
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temporal lattice spacing. We first follow the standard momentum-space construction for 
staggered fermions [22, 26] to identify the eight degrees of freedom that arise from spatial 
doubling (the doubling associated with the time direction is already explicitly taken into 
account in our two-time-slice construction of the transfer matrix).^ We diagonalize the 
Hamiltonian in spin-taste space, and proceed to determine the eigenstates and eigenvalues 
using a generalized Bogoliubov transformation. Two-point correlators can then be easily 
found; they clearly show the expected violations of unitarity in the ghost sector of the 
theory. The explicit calculations below and in App. A demonstrate that Tq has a complete 
set of (right or left) eigenstates in the free theory. In other words, they show that, when 
Uj{x,t) = 1, blocks of the form of Eq. (2.21) do not occur in the Jordan normal form 
of the transfer matrix. As an alternative to the Bogoliubov-transformation approach, a 
path-integral construction of the correlators is given in App. B. 



A. The free Hamiltonian 

Setting Uj{x,t) = 1, writing-^" 

fG4U = l)=exp{-2atH) , (3.1) 

and taking the limit at — > 0, we find the Hamiltonian H in that limit:^^ 

H = 1^2]^ (*i(^)' + Mxf + ni(f)2 + n2(f)') +J2ivj{x) X (3.2) 

X I j 

-<i>i(f)$2(x + /) - <i>i(f + j)$2(x) + n2(x)ni(f + /) + n2(f + /)ni(f) 



Introducing creation and annihilation operators a\, ai, a^, 02 through 


$i(x) = 


y^i=(ai(^) + al(-^))e*% 


ni(f) = 


j^^(a,{k)-a\{-k))e^-^ , 


$2(f) = 


/^H-.)-4(^))e-, 


n2(x) = 


y^^(a2(-fc) + 4(^))e*% 


in which 


f _ f d^k 




Jk~ J (27r)i ' 



(3.3) 



(3.4) 



^ The 16 degrees of freedom are identified as the four tastes of Dirac fermions, each with four spin degrees 

of freedom. 
^° The factor 2 appears because Tq is a double time-shce transfer matrix. 
^^ From now on we drop hats on operators. 
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this can be re-expressed as 

H = -J2 [ [ e'^'^+'">-^ X (3.5) 



2 

m ( a[{-k)ai{i) + ai{k)a[{-£) + al{k)a2i-i) + a2i-k)al{i) 



'k Ji 



+ J2 e*^' e'""^ "^ (-ai{k)a2{-e) + al{-k)al{f) - a2{-k)ai{£) + al{k)a\{-i)') 
j 

Here the factors e ''j are equal to the sign factors ri'j{x) in Eq. (3.2), if we choose 

7?,,/ = (tT, TT, Tt) , TT^/, = (0, TT, Tt) , TT^^ = (0, 0, Tt) . (3.6) 

The creation and annihilation operators have commutation rules 

K{k),al{i)] = 6{k-i)6^p. (3.7) 

We now split up the (spatial) Brillouin zone as in Refs. [22, 26] for staggered fermions:^^ 

k =p + 7iA , i = q + TiB , (3.8) 

with 

v?A,7fBG{(0,0,0),(7r,0,0),...,(7r,7r,7r)}, (3.9) 

such that — 7r/2 < Pj,qj < 7r/2. Operators get relabeled as in 

aa{k) = aa{p + if a) = a'^ip) , (3.10) 

etc. Performing the sum over x in Eq. (3.5) we find the delta functions 

(5(p + g + VTA + vtb) = S{p + q)6AB, (3.11) 

S{p + q + tta + ttb + vf^p = 6{p + q)X^^^ , 

where the second of these equations defines three symmetric matrices X^ . In the Hamiltonian 
these matrices occur in combination with the factors e'^^^^ coming from e*^^, and we define 
another set of matrices 

a\B = X':,^e^--^ . (3.12) 

The matrices a^ are real, and antisymmetric: 

4a = ^A e'-'^ = X^AB <^'^"^^'^' = -X'ab e^^"^ = -<B , (3.13) 

because (vr^')j — ^- Hence the a^ are anti-Hermitian. Using all this, we can simplify the 
expression for the Hamiltonian to 

H = / S 2"^ {ai{p)ai{p) + ai{p)a\{p) + 4(p)«2(p) + a2(p)4(p)) (3-14) 

+i y^ sin (pj) ( ai{p)a^ a2{p) — a2a^a\{p) 



^^ See also Ref. [2] for a review. 
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Note that the integral over p is over the reduced Brillouin zone. 

Finally, the eigenvalues of the 8x8 matrix ^ . sin {pj)a^ are equal to ^isi^p) with s^{p) = 
Y2j sin^ (pj). Dropping a constant proportional to 1/a^, we can thus write if as a sum and 
integral over terms of the form 

h{p) = m (al{p)ai{p) + al{p)a2ip}) ± s{p) (ai{p)a2{p) - 4(p)4(p) ) • (3-15) 



B. Eigenvalues and eigenstates 

Dropping the dependence on p in Eq. (3.15), our next step is to find eigenvalues and left 
and right eigenstates of the non-Hermitian Hamiltonian 

h = m{a\ai + a\a2) + s{aia2 — a2'^i) , (3.16) 

in which Oi 2 and a\ ^ are a set of bosonic annihilation and creation operators, m > and s 
is real. (Taking s real without restriction on its sign takes care of both signs in Eq. (3.15).) 
Adapting the method of Ref. [27], we introduce new operators 

h\ = cos9 ai — sinO al , (3-17) 

62 = cos 9 a2 — sin 9 a\ , 

61 = cos 9 a[ + sin 9 02 , 

62 = cos 9 02 + sin 9 ai , 

which obey the commutation rules 

[baM=SaP^ (3.18) 

while all other commutators vanish. Note that the operators ba are not the Hermitian con- 
jugates of the operators ba, which is why this is a "generalized" Bogoliubov transformation. 
Expressed in terms of these operators h becomes 

h = {mcos{29) + ssin(2^)) (bibi + M2) (3.19) 

+ (— msin(2^) + scos(26')) (6162 — &2&1) + constant . 

Requiring the term proportional to 6162 — ^2^1 to vanish yields 

9 = - arctan(s/r?7,) , (3.20) 

where we picked the solution that vanishes for s — )■ 0. Substituting this solution into 
Eq. (3.19) gives 

h = Vm"^ + s2 {bibi + 62^2) + constant . (3.21) 

Even though ba 7^ &„, the operators A^^ = baba are still number operators, and we can find 
a set of right eigenstates |ni,n2)i? such that 

A^i|'^i,^2)ij = '^i|'^i,^2)i? , (3.22) 

N2\ni,n2)R = n2\ni,n2)R , 
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and 

&i|^i,^2)k oc |ni - 1,^2)/? , (3.23) 

b2\ni,n2)R (X \ni,n2 - 1)r . 

Because of the upper bound (2.20) on the absolute values of the eigenvalues of Tq, which 
implies a lower bound on the real part of the eigenvalues of h, we find that there exists a 
"right vacuum" state |0, 0)/? annihilated by 61 and 62; otherwise ni and n2 could be lowered 
indefinitely, violating this bound. Therefore, it follows that A^^i and N2 vanish on |0,0)ij. 
The right eigenstates of h are then given by 

|n,m)R = ^^6^6-|0,0)r, (3.24) 

h\n, m)fi = (n + m)E\n, m)^ , E = yvrfi + s^ , 

where we have dropped the constant in Eq. (3.21). The normalization of the states we have 
chosen is convenient but arbitrary, since the only relevant normalization condition relates 
right to left states (Eq. (3.26) below). A similar reasoning leads to a left ground state l(0, 0|, 
and a construction of the left eigenstates 

i(n,m| = ^(0,01—^6^6™, (3.25) 

\/n\vn\ 

L{n,m\h = L{n,m\E{n + m) , 
where the normalization here follows from Eq. (3.24) if we demand that 

L{ni,mi\n2,m2)R = 5„i,„2 Smum2 ■ (3-26) 

We note that 

i(n,m|^(|n,m)^)t (3.27) 

because 6q, 7^ fej^,; there is no simple relation between left- and right-eigenstates. However, 
we do have a completeness relation: 

'^\n,m)n L{n,m\ = l . (3.28) 

n,m 

Completeness is not obvious, since our Hamiltonian is not Hermitian. We refer to App. A 
for a proof. 

Under PT symmetry, /i(p, s) -^ h{—p, —s), so it is actually the sum h^p, s) + h{—p, —s) 
that is PT symmetric. Since, from Eq. (3.21), h{p, s) and h{—p, —s) have identical eigenval- 
ues, PT symmetry implies that those eigenvalues must either be real or come in complex- 
conjugate pairs. In fact all the eigenvalues are real, as seen in Eq. (3.24), and all energy 
eigenstates of the sum may be chosen to be eigenstates of PT. In the PT literature [25] 
this situation is referred to as "unbroken PT symmetry." This is somewhat different from 
standard field-theory usage of the terms broken and unbroken symmetry, which refer to the 
properties of the ground state only. 
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C. Two-point correlators 

The free ghost partition function Z is the trace of the (Euchdean) evolution operator, 
which we may write, using Eq. (3.28), as 

Z = J2L{nMe-^'\n,m)n = J^e'^-^-)^^ = _ ^ . (3.29) 

n,m n,m ^ ' 

From Eqs. (3.24), (3.25) and the commutation rules (3.18), it is straightforward to show 
that (with no sum over a) 

{ha{t)Km = -J2L{n,m\e-^^-'Ho,e'%^\n,m)n = ^ _ ^_^^ , (3.30) 



n,m 



and likewise 



1-e 
where we used 



{bait^m = :z^ , (3.31) 



oo 



d sr^ ^c-r d 1 e 



-ET 



ho d{ET)h^' d{ET)l-e-^T ^i_,^ETf ^''■^'^ 

More interesting are the two-point correlation functions involving the original creation and 
annihilation operators a^ and a]^. For these we find from Eqs. (3.30) and (3.31), using the 
inverse of Eq. (3.17), 



Et I p-E{T-t) 



s e + e 



and 



(al(t)a,(0)) = 5., -— — — -^ + — — — — ^ ) . (3.34b) 



a,{t)a,m = <5,+,-3 — ^ _ ^„^^ = -{al{t)a]m , (3.33) 



/ /x t/^Nv r I'E + m e"^* E-m e'^^'^'^^ \ /„„. x 

{a^{t)a]iO)) = 5,, ( ^^ ^-^3^ - ^^ Y^-r^) ' (3-34a) 

E-m e"^* E + m e'^^^^"*) 

2E 1 - e-^^ ^ 2E 1-e- 

Equation (3.34b) is a clear indication of the violation of unitarity in this theory: In the limit 
T — 7- oo, this correlator is negative (for s ^ 0). In a normal theory, it would be a sum of 
decaying exponentials times positive coefficients. 

An alternative, path-integral derivation of the correlators Eqs. (3.33) and (3.34) is given 
in App. B. That approach also makes possible a direct comparison between our treatment 
of this nonunitary theory, and a treatment of a similar Hamiltonian in the PT symmetry 
literature [28], where unitarity is restored through a redefinition of the theory. 

IV. THE EFFECTIVE THEORY 

We now return to the interacting theory defined by the transfer matrix Ttotai given in 
Eq. (2.25). In order to derive the chiral effective theory for PQQCD, we ffist need to 
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identify the light degrees of freedom, which we will collectively refer to as "pions." Here 
an extension of the Vafa-Witten theorem [18] turns out to be useful, as already observed 
in Refs. [17, 29]. This is discussed in Sec. IV A. Further, in the case of unquenched QCD, 
Lorentz invariance plays an important role [16]. Our setting is Euclidean, and we instead 
have hypercubic invariance, which we assume to enlarge to 0(4) in the continuum limit, 
as usual. We use this to argue that our pions satisfy the expected dispersion relation in 
Sec. IV B. 

We then put all the ingredients together to derive the effective theory in Sec. IV C 
A subtlety absent in the unquenched theory arises because, in our theory with arbitrary 
fermionic and ghost quark masses, pion masses can in general be complex. In Sec. IV D, 
we use CPT invariance, already introduced in Sec. II C, in order to show that the phases of 
LECs in our theory are the same as in normal ChPT. This does not preclude the occurrence 
of complex masses in the effective theory, but the effective theory can now be used to 
investigate this issue in more detail, as we do in Sec. IV E. 

As before, we will assume that all quark masses are always positive, both in the mass- 
degenerate and mass-nondegenerate cases considered below, as this is the setting for which 
the path integral (2.16) is well-defined, and the arguments of Ref. [18] apply. For simplicity, 
we will also assume that the continuum limit has been taken, so that we can ignore the 
peculiarities of the staggered quark formalism with respect to species doubling. In particular, 
our arguments will apply to QCD with any number of continuum quarks. 

A. Pion spectrum: extension of the Vafa Witten theorem 

PQQCD contains unquenched QCD. Concretely, this means that if we consider correlation 
functions of operators made out of sea quarks (and gluons) only, these correlation functions 
coincide exactly with those of QCD with only sea quarks, i.e., unquenched QCD. Therefore, 
we know that PQQCD has excitations which correspond to pions made only out of sea 
quarks, and that all correlation functions made out of sea-pion operators behave as they 
should in a healthy quantum field theory. 

First, consider PQQCD in which all quark masses are equal; i.e., rris = m^ = irtg.^^ In this 
case, we have a vector flavor symmetry group^^ SU{Ns + Ny\Ny), which, if it is unbroken, 
relates the two-point (and other) correlation functions of all pions in the theory to each 
other, and thus implies that there is a fully degenerate multiplet of pions in the adjoint 
representation of SU{Ns + N^\Ny). As already observed in Ref. [17], the vector- flavor group 
SU{Ns + A^^|A^^) is unbroken because of an extension of the Vafa-Witten theorem [18]. 

Reference [18] contains two proofs. The first proof is based on a consideration of quark 
bilinears, and goes through without modification in the partially quenched case. It implies 
immediately that valence condensates are equal to sea condensates, because with m„ = rris 
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For the rest of this article we wiU use the subscripts s, v, and g to refer to sea, valence, and ghost quarks, 

respectively. 

Until further notice, we use for simplicity the language of the "fake symmetries" introduced in Ref. [1], 

which do not take into account all the subtleties coming from the ghost sector, rather than the correct 

symmetries introduced in Refs. [17, 19, 21]. The subtleties do not affect the general argument. In Sec. IV D, 

where the details of the chiral Lagrangian matter, we use the nonperturbatively correct form. See Ref. [2] 

for a review of the issues arising from the ghost sector. 
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there is in fact no distinction between these two types of fermions. One ingredient that is 
needed is that the measure in the path integral is positive, but this is not changed by the 
fact that the valence part of the fermion determinant is missing. It is also straightforward 
to prove that no flavor-symmetry breaking can take place between the valence and ghost 
sectors, at the level of order parameters made out of quark bilinears [29]. 

As Ref. [18] points out, spontaneous breaking of flavor symmetry could occur without 
quark-bilinear order parameters. A more general proof that it does not considers the current- 
current correlation functions for the conserved flavor currents. If spontaneous symmetry 
breaking took place, these correlation functions would couple to the corresponding massless 
Goldstone excitations, and would therefore show a power-like fall-off. The idea is to show 
that such correlation functions instead satisfy an exponential bound that is uniform in the 
gauge- field configuration. This implies that no Goldstone mesons exist, and thus that flavor 
symmetry is unbroken. Since this argument is based on the Ward-Takahashi identities for 
flavor symmetry, the framework also applies to the Euclidean partially quenched theory, 
because these identities can also be worked out from the Euclidean path integral. 

In fact, Ref. [18] establishes a somewhat more complicated bound by considering a 
smeared, gauge-invariant quark propagator. But the key point for our discussion is that 
this bound is obtained for a fixed gauge-field background, and independent of that gauge- 
field background. Since the only difference between the standard QCD case and the PQQCD 
case is the relative weight of all gauge-field backgrounds in the Euclidean path integral, the 
analysis of Ref. [18] carries over to the partially quenched case. 

Therefore, we conclude that the Vafa-Witten theorem applies to the PQQCD case with 
fully degenerate quarks [17]. PQQCD contains a complete SU{Ns + N^lN^) multiplet of 
Goldstone mesons, in the mass-degenerate case. As in full QCD, if there are no other 
excitations that are accidentally light, one can consider the low-energy regime, in which 
this partially quenched pion multiplet contains the only light excitations below a certain 
scale. Since we know that the pions in the sea sector have a nonzero mass, all pions in the 
partially quenched theory are massive, and thus all correlation functions of pion operators 
fall off exponentially, with a rate equal to the pion mass. 

In the partially quenched theory, there are other states, made for example from multiple 
valence quarks, that have no analogue in the sea sector, so these states are not constrained 
by the Vafa-Witten argument above. The absence of spontaneous symmetry breaking means 
that there is no fundamental reason for these states to be light, but we cannot eliminate the 
possibility that they are accidentally light or massless. As in Ref. [16], this possibility can 
be excluded only by assumption: we assume that the pions are the only light states in the 
theory. 

With this additional assumption, all correlation functions must decay exponentially, 
which we may call "cluster-like." This result is supported by overwhelming numerical evi- 
dence from lattice QCD computations. For true clustering, it is necessary in addition that a 
"vacuum" state exists and that this state is nondegenerate. The existence of a complete set 
of states follows from the existence of a transfer matrix for the partially quenched theory, 
cf. Sec. II C.^^ A state with a maximal absolute eigenvalue of the transfer matrix must also 
exist, because the transfer matrix is bounded. We cannot, however, prove uniqueness. For 



^^ At least this is true at nonzero lattice spacing; we will ignore subtleties with defining a Hilbert space in 
the continuum limit. We also will follow Sec. II B in assuming that any gauge configuration on which Tq 
does not have a complete set of eigenstates can be ignored. 
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example, we must by assumption exclude the possibility that some breaking of a discrete 
symmetry {e.g. parity) occurs, resulting in two vacua, but without Goldstone bosons. With 
these assumptions, which are closely analogous to the assumptions required in Ref. [16], it 
follows that the degenerate partially quenched theory obeys clustering. 

Next, we consider the nondegenerate case, always keeping all quark masses nonzero. 
When we move away from the degenerate point by taking the valence, ghost, and sea masses 
unequal (including, but not limited to, the partially quenched case where valence and ghost 
masses remain degenerate), we can still apply the bounds of Ref. [18] on correlation functions 
directly, even if the vector symmetries are broken explicitly by the mass differences. The 
gauge measure remains positive in this case because it is just a product of two positive 
determinants (valence and sea determinants) divided by the positive ghost determinant. 
Thus the bounds go through, and all connected correlation functions from point x to point 
y decay exponentially. Here "connected" means they are formed out of one or more quark 
propagators that go from x to y. For connected correlation functions, we thus automatically 
have behavior which is "cluster-like," in that correlators decay exponentially. 

However, the partially quenched theory has many disconnected correlators, most of which 
have no pure-sea analogue (made using many valence flavors so that x to y contractions 
do not occur). We again need to assume that such correlators do not have power-law or 
anomalously light decay. With this assumption, the nondegenerate partially quenched theory 
has exponential decay in all channels. If we further assume, as above, that there is a unique 
lowest state, then the partially quenched theory obeys clustering, even with nondegenerate 
masses. 



B. The dispersion relation 

We will assume that the continuum limit of the Euclidean partially quenched theory has 
0(4) "space-time" invariance. What we wish to argue next is that this leads to the expected 
form of propagators for one-particle pion states, and thus to the usual relation between 
energy, mass and spatial momentum. 

The transfer matrix for the full theory has (right) eigenstates, which can be classifled 
according to their eigenvalues. In addition, since the theory is invariant under spatial trans- 
lations, spatial momentum is conserved, and the eigenstates of the transfer matrix can 
simultaneously be labeled by their spatial momentum. This follows because the transfer 
matrix, even if it is not Hermitian, generates translations in the time direction, and thus 
commutes with the generators of spatial translations, i.e., spatial momentum. 

Consider flrst a pion two-point correlator with zero spatial momentum. We know from 
the preceding section that this correlator falls off exponentially for large times: 

C^{t) oc e-"^-l*l . (4.1) 

Here the parameter m^r might in principle be complex, since the Hamiltonian is not Hermi- 
tian, and it is possible at this stage that r/ivr is one of a complex pair of eigenvalues. However 
we do know that 771^^ has a positive real part, as required by the exponential damping of 
correlation functions. There may also be other states that contribute to €-,^{1), which would 
lead to other exponentially damped contributions, with a faster decay rate (for instance, 
three-pion states). 
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The Fourier transform of CT^(t) is 

dtev.t-mAt\^^^^^_ (4.2) 

PI + "ml 

We now consider the correlator of a pion with a nonzero spatial momentum p. By 0(4) 
invariance, the Fourier transform of this correlator has to be equal to 

/(P) = 4^ ' (4-3) 

where p^ = 'YlinPtJi.Pfj'- ^^ ^^ ^'-'^ Fourier transform back, we obtain the leading exponential 
of the pion correlator with nonzero spatial momentum, 



1 [^ , 2m^e-*P** 



m^ 



^'''■'^— J j^^Aw^r^ '''''' <"' 



with E = \/vfi\ +jp. It follows that the dispersion relation for pions in PQQCD is the 
usual one, even though this theory is only defined in Euclidean space and even though the 
"mass" rra^ may in principle be complex. 

C. The chiral effective theory 

We now have all the needed ingredients to extend the arguments of Ref. [16] to the 
partially quenched case. The existence of the transfer matrix, together with the arguments 
and assumptions discussed in Sec. IV A, tell us that the theory clusters, and that the lowest 
states are pions. From Sec. IV B, the pion states have arbitrary momentum, with the usual 
dispersion relation for energy in terms of spatial momentum and (possibly complex) mass. 
Further, from Sec. II, the transfer matrix has a complete set of eigenstates, under the mild 
assumption discussed at the end of Sec. II C. 

We then just follow Ref. [16] step by step. The low-energy theory is dominated by 
exchange of the lightest particles, the pions {i.e., one has pion "pole dominance," although 
of course there are no poles in Euclidean space). Pion interactions in a small region of space 
are described by vertices, which do not depend on the overall process due to clustering. As in 
the normal QCD case, we need to assume that these vertices are strictly local, i.e., that the 
vertices may be expanded in powers of momenta. Long-distance correlation function then 
involve the exchange of pions between local vertices. Since the vertices are strictly local 
building blocks in correlation functions in the effective theory, two (or more) vertices can be 
joined by more than one pion propagator, leading to loops. And because the set of eigenstates 
is complete, the correlation functions can be written as sums over intermediate states, which 
implies that when there are loops we sum over the pion momenta. Here a subtlety could arise 
if the possibility of complex masses is realized, because the standard arguments to relate 
sums over intermediate (on-shell) states to four-dimensional loop integrals would seem to 
depend on analyticity, and the poles and cuts in Green's functions will not be in the normal 
places if there are complex masses. We will assume that no complex masses arise, and then 
note that this assumption is confirmed a posteriori by the effective theory of PQChPT in the 
conventional (and most important) case of degenerate valence-quark and ghost masses. As 
discussed in Sec. IV E, however, complex masses do seem to be possible in PQChPT when 
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valence- and ghost-quark masses are nondegenerate. Thus the reader should keep in mind 
that the foundations of the effective theory are less secure at present if such nondegenerate 
masses are allowed. 

The ingredients above, which we may roughly summarize as "pions, vertices, and loops," 
are all that Ref. [16] needs to argue that there is an effective chiral theory describing the low 
energy behavior of the theory. To get the standard chiral theory (ordinary ChPT in the full 
QCD case, PQChPT in the partially quenched case), a further technical argument is needed 
to show that one may choose the chiral theory to have local chiral symmetry (with vector 
sources inserted that transform like gauge fields, and up to anomalies, which must be dealt 
with separately). That part of the argument is based on the Ward-Takahashi identities. 
Since, from the point of view of Ward-Takahashi identities, the partially quenched theory is 
just like the ordinary theory, but with the chiral group extended to a graded chiral group, 
this part of the argument goes through unchanged. 

D. The effective theory and constraints from CPT 

We work with the nonperturbatively-correct partially quenched chiral Lagrangian intro- 
duced in Refs. [17, 19, 21]. The Lagrangian is a function of the chiral fields S(a;) and T,~^{x). 
S is parameterized as 

EW = exp(2*W//) ; *W = ( ^^^^j' ^g ) (4.5) 

where / is the pion decay constant, and where (f){x) are the quark-antiquark mesons, are 
the ghost-antighost mesons, and u and u are quark- ant ighost and ghost-antiquark mesons, 
respectively. While (p and (p are commuting fields, u and u are Grassman valued. At leading 
order (and in the continuum limit), the Euclidean effective Lagrangian is [19] 

£^g = i^ str (d^J^df,^-^) - V str (MS + S^^ikf^) (4.6) 

8 

where str denotes the supertrace, M is the quark and ghost mass matrix, and f is a LEG. 
We have assumed that the super-?]' field $o = —"i str In E = tr (0 + i(j)) has been integrated 
out, which is possible as long as the theory is not completely quenched [17]. 

In the past, it has been assumed that the LEGs in the partially quenched chiral Lagrangian 
are real, just as the corresponding ones are in the ordinary Hermitian chiral Lagrangian for 
full, unitary QGD. For most LEGs, such as v above, this follows from the fact that QGD is a 
special case of PQQGD, and LEGs that the two effective theories share (for the same number 
of sea-quark flavors) are equal [3]. However there are additional LEGs that are unique to the 
partially quenched theory, and vanish for pure sea quantities [30] . We would like to be able 
to argue that those LEGs are real too, in order to see that complex meson masses will in 
general not occur. ^^ The antiunitary GPT symmetry of the theory, introduced in Sec. 11 G, 
can be used to show this. 

We must first determine the transformation properties of (p, (p, u and u under GPT 
symmetry. This is a bit subtle since the discrete symmetries of staggered quarks (or ghosts) 



^^ They may still occur in the flavor-diagonal sector, as we will see in the next section. 
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include additional discrete taste transformations (see for instance the discussion of spatial 
inversion symmetry in Ref. [26]). Since we want the transformation laws of pseudoscalar 
mesons under continuum CPT (without additional taste transformations), we must look on 
the lattice at taste-singlet mesons. 

From Ref. [31]/^ a staggered quark-antiquark bilinear with spin F and taste S is 

^ ^ g(x + A) q{x + B)ti {n\A) F n{B) Ft) (4.7) 
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A,B 



where A and B hypercube vectors (with components 1 or in each of the four directions), 
F= = S*, VL{x) = 7^^ 72 ^73 ^74'', and we have omitted gauge links for simplicity. For the 
quark-antiquark meson field (pjk, with j and k the quark and antiquark flavors, respectively, 
we have: 

'^^■'^(^) = ^5Z9fe(a: + ^)?,(a; + 5)tr(fit(^)^^Q(5)) 

A,B 
= ^ E (^^+^.^ (-1)^^+^^ q,{x + A) q,{x + B)) (4.8) 

A,B 

where = should be read as "has the same renormalized matrix elements as," D = (1, 1, 1, 1) 
is the diagonal of a hypercube, and y = x + D/2 is defined for convenience to be the middle 
of the hypercube. The taste-singlet pseudoscalar meson is thus created by a "four-link" 
operator that joins opposite corners of a hypercube. The overall factor of i in Eq. (4.8) is 
crucial; it is required to make the propagator {(pjk{y) 4>kj{y')) positive, which follows at the 
chiral level from the definitions Eqs. (4.5) and (4.6).^^ 
For ghost-ghost mesons, the corresponding relation is 

My) = 4 E (^^+^'« (-1)^^+^=^ xl(^ + A) x,(x + B)) . (4.9) 

A,B 

Here there is no factor of i because there is no minus sign in the propagator due to statistics; 
by the definitions in Eqs. (4.5) and (4.6), and have equal propagators. Similarly the 
quark-antiqhost and ghost-antighost mesons are given by 



A,B 

My) = 4 E (^^+^'^ (-1)^^+^^ qk{x + A) xj{x + B)) . (4.10) 

A,B 

^^ The following argument can also be given in the language of Ref. [32]. 

^® The positivity of this propagator in the continuum limit can be most easily seen by noting that the 
corresponding propagator for taste ^5 (taste-pseudoscalar) pions is a sum of absolute squares, which in 
turn follows from the fact that the staggered Dirac operator obeys 2?^ = eDe, where e is a diagonal matrix 
in position space with e(a;) = (_i)^i+^2+a;3+a:4 along the diagonal. The two factors of i in the propagator 
cancel the minus sign from Fermi statistics. In the continuum limit, the propagators for (flavor-charged) 
mesons of arbitrary taste will be equal. 
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CPT 


0(y, 2/4) 


-^ (-^,1/4) , 


CPT 


0(y, 2/4) 


^ -4>^{-y,y4 


CPT 


w(y, 2/4) 


-^ -u;'^(-y,y4 


CPT 


^(f,l/4) 


^ W^(-f,Z/4) , 



Again, no factor of i is needed here to make the {ujjkiy) uokjijl')) propagator have the same 
sign as the corresponding propagator, which is required by Eqs. (4.5) and (4.6). 

We can now determine the CPT transformation rules for the meson fields by using the 
ghost and quark transformation rules, Eqs. (2.33) and (2.35). We find 



(4.11) 



with T the transpose, which acts on flavor indices. The factors of (—1)^'' in Eqs. (2.33) and 
(2.35), combined with the fact that the spin-1/2 fields in the taste-singlet mesons are on 
different time-slices, are important to getting the above signs under CPT. The fact that the 
quark- ant iquark meson field is even under CPT is standard. 
From the definition of $ in Eq. (4.5) we then have 

CPT: ^{x,Xi) ^ -^\-x,Xi) (4.12) 

where t denotes a graded transpose defined on a matrix with commuting diagonal blocks 
and anticommuting off-diagonal blocks by 



C d - I -feT rfT 



(4.13) 



The operation t has the useful property that 

str (E* F* GV . . ir*) = str ([K...GF Ef) (4.14) 

where E,F,G,K are matrices of the form in Eq. (4.13). The property can be proved by 
writing both sides as index sums, counting the number of interchanges of anticommuting 
numbers between the left and right sides, and showing that the sign resulting from the 
interchanges is the same as the sign coming from the factors of —1 in the definition of t. 
Note that, unlike the normal rule for the transpose of a product of matrices with commuting 
entries, this equality is only true under the supertrace. It is not true for the matrices 
themselves, as is already obvious from considering a single matrix instead of a product. 

We can now investigate the consequences of CPT for the chiral theory. Our definitions 
have ensured that antiunitary CPT symmetry leaves the Euclidean action of the fundamental 
theory invariant, so it will also leave the Euclidean chiral action invariant. We will see that 
this implies that the phases of all LECs in the partially quenched chiral Lagrangian are 
the same as they would be in normal ChPT. For a term that is special to the partially 
quenched case (because in the normal case with the given number of sea-quark flavors it is 
not independent due to Cayley-Hamilton relations), the phase of the LEC is the same as it 
would be in the ChPT theory for a QCD-like theory with sufficient numbers of flavors to 
make the term independent. The argument goes as follows: 

1. Since Euclidean rotational invariance requires derivative operators to be contracted as 
dfj,dfj, we may redefine x — )■ —x on the right-hand-side of Eq. (4.12) without changing 
sign of any term. (We postpone discussion of anomaly terms, which may have the 
form e^u\ad^dyd\d„, until item 5, below.) 
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2. Chiral symmetry demands that any term in the Lagrangian be formed from one or more 
supertraces. Therefore Eq. (4.12) and the rule Eq. (4.14) mean that CPT effectively 
interchanges S and E~^ and inverts the order of products of arbitrary numbers S and 
S^^ matrices. The overall transpose on the right-hand side of Eq. (4.14) of course has 
no effect inside a supertrace. 

3. Because CPT is antiunitary it complex-conjugates all LECs. 

4. Thus the effect of CPT on the Lagrangian is exactly the same as the effect of taking 
the Hermitian conjugate would be in a theory where S is unitary. The phases are 
therefore the same as in normal ChPT. 

5. For anomaly terms with e^yXad^dudxd^, an extra minus sign results from taking x — )■ 
—X, since we are not flipping the sign of X4 in Euclidean space. Therefore such terms 
require an overall factor of i (relative to what they would have in Minkowski space) 
in order to be CPT invariant. However, this factor of i is precisely the factor that 
automatically arises in going from Minkowski space to Euclidean space in any term 
with a single time derivative. Thus the LEC again has the same phase as it would be 
ordinary ChPT. 

A consequence of this result is that all parameters in the chiral effective Lagrangian 
corresponding to the square of a meson mass are real (in the usual convention in which they 
would also be real in the case of a normal theory). Further, all mass terms shared by the PQ 
and normal theory must have signs that ensure that squared meson masses are positive. This 
means that as one turns on partial quenching by moving away from a point where valence 
and ghost masses are degenerate with sea quark masses, the squared meson mass terms must 
stay positive, at least until the valence, ghost and sea mass differences become large. ^^ As 
we will see in the next section, this still does not preclude complex masses from appearing in 
the theory for small perturbations, due to subtleties in the flavor-diagonal sector. However, 
"trivial" complex masses, coming directly from complex LECs, are ruled out. 



E. Masses in the chiral theory 

As mentioned in Sec. IV C, the possibility of complex masses is a concern for our deriva- 
tion. Once we have the chiral effective theory, however, we can study this possibility in more 
detail. 

Consider, for instance, the two-point function of a flavor-diagonal valence pion and a 
flavor-diagonal ghost pion, in the mass nondegenerate case. Extending the results of Ref. [1], 
this two-point function is given by 



>,(p2 + M2)(M2 - M2) + Ns{p^ + M2)(p2 + M2) 



19 Negative squared masses would not imply imaginary masses, but rather signal a phase transition and the 
need to find a new vacuum state. We believe such transitions, driven by possible mass terms special to 
the partially quenched theory, cannot occur in the usual PQ case of degenerate valence and ghost masses, 
since the dynamics of such a theory are controlled by the sea masses. We have not however ruled out the 
existence of such transitions when valence and ghost mass differences are large. 
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where we have chosen the number of ghost flavors, A''^, equal to the number of valence 
flavors, Ny, while allowing rrig ^ my, and where we have sent the singlet part of the rj' mass 
to infinity. M^^g^s are the masses of (flavor nondiagonal) valence, ghost, and sea pions; in 
this subsection we use capital letters to distinguish meson masses from quark masses. The 
denominator of this expression is a quadratic form in —p"^ of the form 

Ai-pY + Bi-p^) + C , (4.16) 

with 

A = Ns, (4.17) 

B = -Ny{Ml - Ml) - Ns{Ml + Ml) , 
C = NyMliM^ - Ml) + NMlMl . 
The discriminant of this quadratic form is 

B^ - AAC = {Ny + Ns)\M^ - Mlf + 4Ar,Ar,(M| - MI){MI - Mf) . (4.18) 

This vanishes for Mg = My, and Eq. (4.15) has a double pole in that case, as we expect, 
unless also Mg = My. 

In order to obtain only real single poles, B^ — AAC needs to be positive, and this turns 
out not to be always the case. For simplicity, take Ny = Ns = N, so that 

{B^ - AAC) \n.=n.=n = AN\Ml - Ml){M^g - M^ . (4.19) 

This is only positive if Mg < My and Mg < M^, or if Mg > My and Mg > M^. This means 
that even a small perturbation from the degenerate case can lead to a situation in which the 
discriminant is negative, leading to a conjugate pair of complex zeros of Eq. (4.16). The real 
part of these zeros is given by —B/{2A) = Mg. The two-point function is still well defined, 
but the effective theory implies that it is possible to be in the situation that energies of 
states become complex (with a positive real part). Note that the existence of a conjugate 
pair of energies is consistent with the PT symmetry of the theory. The reality of Mg, Ml 
and Ml, which follows from the reality of the LECs of theory, does however not guarantee 
that the poles of this real propagator occur at real values of —p^. 

Another instructive example is to take My = Ms 7^ Mg (while leaving Ny and A^s arbi- 
trary), which keeps B^ — AAC > 0. Now the zeros — p^ of the quadratic form are 

-p'_ = Ml, (4.20) 

-pI = M^ + fiM^g - Ml) . 



Ml>M^(l + -^)> Ml . (4.21) 



However, the zero —p\ becomes negative for 

Ny 

This would make the Fourier transform of Eq. (4.15) ill-defined, but this can only happen if 
one perturbs My far enough away from Mg. This suggests that for a choice of My such that 
Ml = M^(l + Ng/Ny) a phase transition takes place. In that case, the effective partition 
function would have to be evaluated by performing a different saddle-point expansion than 
the one assumed in writing down Eq. (4.15). Of course, as already remarked above, no such 
problems, and no complex energies, appear when Mg = My, which is the "physical" case of 
PQQCD, because in that case 5^ — AAC = 0, and we recover the usual double pole. We 
therefore do not pursue the possibilities arising for Mg 7^ My further. 
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V. CONCLUSION 

In this article, we have presented the theoretical evidence that PQChPT provides the 
correct low-energy effective theory for PQQCD. Our starting point is the discussion by 
Leutwyler [16] of the foundations of chiral perturbation theory in the case of full QCD. The 
cluster property of a Lorentz invariant, local quantum field theory plays a central role in 
Ref. [16]; unitarity, in contrast, appears not to be needed. This starting point is essential in 
any attempt to extend the validity of ChPT to the partially quenched case, in which unitarity 
is lost. Therefore, our main task has been to see to what extent the cluster property can be 
established in PQQCD as well. 

The key ingredients are the existence of a bounded transfer matrix, as well as an extension 
of the Vafa-Witten theorem to the partially quenched case. The existence of the transfer 
matrix allows us to identify a complete set of states in the theory, even though the transfer 
matrix and the corresponding states do not have all the usual properties they possess in a 
unitary theory. An important point, however, is that the pion states do satisfy a rotationally 
invariant dispersion relation. The Vafa-Witten theorem [18] allows us to connect the space- 
time dependence of correlation functions in the partially quenched theory with those of 
correlation functions of the corresponding full QCD theory with the same set of sea quarks. 
While these ingredients together are not sufficient to provide a rigorous proof of clustering, we 
have identified the further assumptions needed to establish the cluster property in PQQCD, 
assuming that it holds in the full theory. ^*^ 

Once the cluster property is established, the argument for the correctness of PQChPT 
as the low-energy effective theory for PQQCD follows mainly along the same lines as that 
given in Ref. [16]. While in both cases some additional assumptions are needed (such as 
the assumption that no other accidentally light states exist in the theory), these additional 
assumptions in general have little to do with the "sickness" of the partially quenched theory, 
and we thus believe them to be equally plausible in both the full and partially quenched 
cases. One issue unique to the partially quenched case is the possibility of complex energies, 
which could make it difficult to write effective-theory loops as normal four-dimensional 
integrals, and require instead three-dimensional integrals over on-shell states. The arguments 
presented in Sec. IV support the assumption that this situation does not occur in the only 
case of practical interest, in which ghost quark masses are chosen equal to the corresponding 
valence quark masses, i.e., the case usually referred to as partial quenching. Section IV D 
shows that CPT symmetry requires the effective theory to have real low energy constants. 
This restricts the problem of complex masses to the flavor-diagonal sector. As demonstrated 
in Sec. IV E, complex masses (or poles located at negative values of (Euchdean) — p^) do 
not occur when all ghost quark masses are chosen equal to the corresponding valence quark 
masses. The only "sickness" is the familiar occurrence of double poles, at positive values of 
— p^, in the valence sector. 

Our framework for the construction of the effective low-energy theory generalizes to the 
fully nondegenerate case, in which ghost-quark masses are not equal to valence-quark masses. 
It turns out that this generalized theory contains some new properties absent in the par- 
tially quenched case; these properties can be investigated in the corresponding effective field 
theory. In particular, one can choose values for ghost and valence masses such that complex 



^° We are not aware of a rigorous proof of the cluster property for full QCD either; for a discussion of the 
cluster property in quantum field theories, see Ref. [33]. 
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poles appear in the disconnected propagator. We leave open the question of the validity 
of the standard loop expansion of the effective theory for this case. Of course, as already 
emphasized above, this issue is primarily of academic interest, since, by construction, ghost 
and valence quark masses are always equal in numerical applications of PQQCD. 
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Appendix A: Completeness 

In this appendix, we prove completeness in the free theory of the basis of right eigenstates 
given by Eq. (3.24), or the left eigenstates given by Eq. (3.25). The argument has three steps. 
First, we prove that |0,0)r, l(0,0| are the unique right and left vacuum states. We then 
show that all right eigenstates are obtained by acting with the raising operators 6i and 62 
on |0,0)ii; (and similarly for left eigenstates). Finally, we prove that that no blocks of the 
form of the form (2.21), or larger generalizations thereof, occur in the Jordan normal form 
of the Hamiltonian matrix of the free theory. 

The uniqueness of the right and left vacua can be proved by working in position space 
and solving the differential equations corresponding to the definitions of |0,0)ij and i(0,0|, 
which are given in operator form by 

bi\0,0)R = (cos e ai- sin 6 4) |0, 0)ij = , (Al) 

62|0,0)r = (cos0 aa-sin^ a\)\0,0)R = , 

and 

^(0,0|6i = l(0,0| (cos^ al + sin^ 02) =0 , (A2) 

^(0,0|62 = l(0,0| (cos^ 4 + sin^ ai) = , 

where we used Eq. (3.17). 

Using X for the position associated with oscillator described by ai, a[ and y for the 
oscillator described by 02, al, we have 

(A3) 



1 / 9 ^ 


^'^=72^"" 


d 

dx 


1 , d^ 


"^=^2^"" 


d 
dy 



Now let 

{x,y\0,0)R = ^{x,y), l(0, 0|x,y)/j = x*(x,?/) (A4) 
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be the position-space wavefunctions for the vacua. The equations for ^ip^x, y) and x(a;, y) are 

d 



cos 9 [x + 
— sin 9 I X ~ 
cos 9 i X + 
sin 9 I X — 
The solutions are (recalhng Eq. (3.20)) 

ip{x^y) = Aexp 

x{x,y) 



d_ 
dx 

d_ 
dx 

d_ 
dx 

d_ 
dx 



— sin 9 \y — 

-\- cos 9 [y + 
+ sin 9 [y ~ 
+ cos 9 [y + 



E 



dy 
d_ 

dy 
d_ 

dy 
d_ 

dy 



i'i.x.y) 
^{x,y) 

x{x,y) 
x{x,y) 



0, 
0, 
0, 
. 



(A5) 



-— (x^ + y^) H xy 

2m m 



Aexp ( ^^(^^ + y'^) ^ :^^y 



(A6) 



Since the only free parameter in these solutions is the normalization A, we have shown that 
the right and left vacua are unique. 

It is now clear (following the argument after Eq. (3.23)) that any right eigenstate of h can 
be repeatedly lowered by application of hi and 62 until we reach |0, 0)ij. This is all we need 
in order to show that the eigenstate can, in turn, be obtained by operating with hi and 62 
on the vacuum, and is therefore just proportional to one of the states |n,m)ij in Eq. (3.24). 
For example, consider an eigenstate |l',0')ij, with 



6i6i|l',0% = |1',0% 
M2|l',0% = 0. 



(A7) 



By uniqueness of the vacuum. 



6i|1',0%(x|0,0)r. (A8) 

Operating on both sides of this equation with hi and using Eq. (A7) then proves that 

|l',0')i?(x6i|0,0)^, = |l,0)R. (A9) 

It is straightforward to extend this step to a proof by induction that, similarly, 

|n',m% oc \n,m)R , (AlO) 

for an eigenstate \n',m')fi of the number operators hihi and 62^2 with eigenvalues n' and m' 
respectively. An analogous argument extends this first step to left eigenstates as well. 

The fact that we have shown that all right and left eigenstates are of the form of 
Eq. (3.24) and (3.25) does not yet quite prove completeness of these states. Since the 
Hamiltonian (3.21) is not Hermitian, it is possible that its Jordan normal form has non- 
trivial blocks of the form (2.21), with k ^ 0. We now prove that this does not happen in 
the free theory. First, the Hamiltonian is proportional to the sum of two number operators. 
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A^i = bibi and N2 = &2&2j which commute with each other. Therefore, the matrix repre- 
sentation of the sum is a direct product of the two matrices representing A^i and A^2- It is 
thus sufficient to show that no nontrivial Jordan blocks occur in the matrix representation 
of iV = bb, where b = bi and b = bi, with i = 1 or 2. 

Now assume that a hypothetical Jordan block of arbitrary size p x p occurs, with p > 
1. Such blocks have an eigenvalue on the diagonal, and an arbitrary constant k, on the 
"superdiagonal." The block is trivial ii k = 0; nontrivial otherwise. We already know that 
possible eigenvalues of the number operator N are n = 0, 1, 2, . . . . So, for example, a block 
for eigenvalue n and p = 4 would look like 



/ n K 0\ 

n K, 

n K, 

\0 n J 



(All) 



We will prove by induction on n that no such blocks can exist for any n, and any size p. For 
any such hypothetical block, let \n)ji be the true right eigenvector and |A*^"))^ be the ffist 
nontrivial generalized right eigenvector. These states obey 



N\n)ii = n\n)ii , 
For example, in the p = A case above, \n)ji and \A'^"-^)r could be represented by 



(A12) 



/1\ 






respectively, with 7 arbitrary. 

The theorem that we will prove is: 



and 



/7\ 

1 





(A13) 



If, for any n, there exists a state 1^*^"^)^ that is linearly independent of 
the eigenvector In)^ and obeys 



iV|A("))^ = n|A("))R + «:|n)fi 



(A14) 



then K, = 0. 



Thus \A^^'')ji is a right eigenstate degenerate with \n)ji. However, since we showed above 
that the eigenstates of A^ are nondegenerate, \A^"'^)ji would have to be proportional to \n)ji, 
which contradicts the assumption of linear independence. Thus, no state \A^"'^)ji can exist, 
and the Jordan block is not only trivial {k = 0) but is actually one-dimensional {p = 1). 

We ffist show the theorem is true for n = 0. For a nontrivial Jordan block with n = 
we would have 

N\A^'^^)r = k\0)r, (A15) 



which implies 



Nb |A(°))^ = {bN - 6)|A(°))^ = -b |A(°)) 



R 



(A16) 
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So b\A^^^)ji would have to be an eigenstate of A^ with eigenvalue —1, which is impossible 
since A^ is positive semidefinite. Thus b \A^^^)r must vanish, which implies 

iV|AW)R = 66|A(°))« = 0. (A17) 

Comparing with Eq. (A15), we see that k = 0, as desired. 

We now assume that the theorem is true for eigenvalue n, and prove it is also true for 
eigenvalue n + 1. Suppose we have a state \A^'^~^^^)r that is linearly independent of \n + 1)r 
and satisfies 

N |A("+i))« = {n + l)\A'~^+'^)n + «:|n + 1)r . (A18) 

Then 

Nb I A("+i))^ = {bN - 6) I A("+i))^ = nb I A("+i))^ + V^^^Tl i^\n)R. (A19) 

So the state b\A^"-~^^^) ji satisfies Eq. (A 14) for eigenvalue n, with k replaced by k' = y/n + 1 k. 
Furthermore, b\A^"'~^^'>)ji is linearly independent of \n)ji. If not, b \A^'^~^^^)]i = a\n)ji for some 
constant a, which implies by Eq. (A18) that 

bb |A("+^))r = {n + 1)|A("+^))r + K\n + 1)r = ba\n)R = V^^+l a\n + 1)r , (A20) 

contradicting the assumption that \A^'^'^^^)r is linearly independent of \n + 1)r. Therefore, 
b \A^'^'^^^)ji satisfies the conditions of the induction hypothesis, which implies k' = 0, and 
hence k = 0. 

This concludes the proof of our theorem, and hence completeness of the free theory as in 
Eq. (3.28). 

Appendix B: Path integral for the free ghost theory 

Starting from Eq. (3.16), we reconstruct a path integral which yields the same correla- 
tion functions as those generated by the transfer matrix defined in terms of h. This is an 
alternative to the operator analysis based on the generalized Bogoliubov transformation, 
Eq. (3.17). For simplicity, we take T — > cxd from the beginning. 

We start by defining new fields ipi and conjugate momenta pi 



ijji = —=={ai + a]) , pi = -i^ —{ai-a]) , (Bl) 

yzm V ^ 

which are conventionally normalized, 

bPt,Pj] =iSij . (B2) 

The Hamiltonian then may be written 

h = - {pI + pI + m'^{ipl + ^D) + is {^iP2 + ip2pi) + constant . (B3) 

Introducing sources Jj for ipi and Ki for pj, the partition function in the Hamiltonian form 
of the path integral is 

Z[Ji, Ki] = / Vipi Vip2 'Dpi Vp2 exp{ipj'tpj - h + Jjipj + KjPj) , (B4) 
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with a sum over the repeated index j, and with ipj = dipj/dt. The quadratic integrals over 
pi and p2 are easily done,^^ and the result is 



Z[Ji,Ki\ = jvi,^Vi,2e 



CE[J^,Ki_ 



(B5) 



with the Euclidean Lagrangian £^[Jj,i^j] given by 



CE[Ji,Ki 



^2 ^ ^2 ^ ^2^^2 ^ ^2) _ ^^^(^^ _ ^^^) _ ^^^^^^ _ ^^^) 



-JlV'i - J2i^2 



\{Kl + K, 



sj^i^M 



(B6) 



The terms quadratic in Kj are standard, and just give contact terms that we drop from 
now on. We will also drop the final, total-derivative term, since we want to calculate 
a partition function with periodic boundary conditions on ipj. Note, however, that for 
computing transition amplitudes between tpj eigenstates, dropping the total derivative would 
be incorrect; indeed that term leads to violations of unitarity. Violations arise because the 
Minkowski-space transition amplitude 



Tba = {i^l,b'ip2,b;tb\lpl,a'ip2,a;ta) 



(B7) 



is proportional to exp[s{ipi^b4'2,b — 'ipi,a4'2,a)] and does therefore does not obey Tba = T*^. 

If we focused only on correlation functions of ipj, the sources Kj could also be dropped 
completely, and one would immediately see that the correlation functions are the completely 
normal correlators of two harmonic oscillators with frequency E. One could in fact change 
the rules at this point, and write down a standard (double) harmonic oscillator Hamiltonian 
h' that gives these same ifjj correlators: 



h'='^{pl + pl + E\^^l + ^l)) 



(B8) 



The analogous step for the Hamiltonian of Ref. [27] is taken in Ref. [28]. However, this 
is really a change of rules, because it changes the correlators involving pj. As we will see 
below, nonunitary effects we already obtained in Sec. HI C show up in those correlators. 

From Eqs. (B5) and (B6), we can easily compute all two-point correlators by taking 
derivatives with respect to Jj and Kj. The results are: 



{ut)pm) 

{P^{t)pM) 



-Et 



Sij 
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-%se 



-Et 



^i+j,3 



—le 



-Et 
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—ise 



-Et 



i+i,3 



2E 



rri^e ^* 



2E 



-Et 



5i 



2E 






—le 



Et 



-Et' 



2 



(B9a) 
(B9b) 
(B9c) 
(B9d) 



The Pj have open boundary condition, as is always the case in the path integral formulation. 
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Here, the corresponding results in the standard theory defined by h' are shown in square 
brackets. For both our theory and the standard theory, these correlators satisfy the condi- 
tions required by the commutators Eq. (B2) in the limit t — )■ 0. (The relevant correlators 
are unaffected by the dropped contact terms.) 

Aside from the off-diagonal terms in the ip-p correlators, the difference between our 
theory and the standard theory appears in the normalization of the p-p correlator in 
Eq. (B9d). From Eq. (Bl) and Eq. (B9), we can easily check that we reproduce the correla- 
tors Eqs. (3.33) and (3.34) for the creation and annihilation operators in the limit T — >■ oo. 
In particular, the unitarity- violating negative result for (aJ(t)aj(O)) shows up here because 
of the incomplete cancellation of the various contributions from Eqs. (B9a) through (B9d). 
(Note that, in order to recover expected results for correlators of creation and annihila- 
tion operators in the standard theory, one must replace the definitions Eq. (Bl) by letting 
m — )■ £", so that the Hamiltonian h' takes its standard form in terms of a^ and a\.) 
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